% function [baseMVA, bus, gen, branch, genpara, NTswitch, sbspara,Hess_coef,u0,mu_coef] = case39
baseMVA = 100.00;

% bus  type1 type2  Pdemand  Qdemand   Gshunt   Bshunt   Vmodule    Vangle  baseKV  Vmax  Vmin
% type1: PQ bus = 1; PVbus = 2; reference bus = 3; isolated bus = 4
% type2: non-load = -1; load -- no shedding = 0; load -- shedding =1
bus = [
     1     1    -1      0.0000      0.0000      0.0000      0.0000      1.0540     -7.6511      0.0000      1.0600      0.9400
     2     1    -1      0.0000      0.0000      0.0000      0.0000      1.0448     -4.2809      0.0000      1.0600      0.9400
     3     1     1    322.0000      2.4000      0.0000      0.0000      1.0417     -7.8674      0.0000      1.0600      0.9400
     4     1     1    500.0000    184.0000      0.0000      0.0000      1.0368     -9.3334      0.0000      1.0600      0.9400
     5     1    -1      0.0000      0.0000      0.0000      0.0000      1.0490     -8.4400      0.0000      1.0600      0.9400
     6     1    -1      0.0000      0.0000      0.0000      0.0000      1.0537     -7.8430      0.0000      1.0600      0.9400
     7     1     1    233.8000     84.0000      0.0000      0.0000      1.0410     -9.8256      0.0000      1.0600      0.9400
     8     1     1    522.0000    176.6000      0.0000      0.0000      1.0388    -10.2709      0.0000      1.0600      0.9400
     9     1    -1      0.0000      0.0000      0.0000      0.0000      1.0541     -9.9292      0.0000      1.0600      0.9400
    10     1    -1      0.0000      0.0000      0.0000      0.0000      1.0600     -5.8866      0.0000      1.0600      0.9400
    11     1    -1      0.0000      0.0000      0.0000      0.0000      1.0567     -6.5550      0.0000      1.0600      0.9400
    12     1     1      8.5000     88.0000      0.0000      0.0000      1.0434     -6.5634      0.0000      1.0600      0.9400
    13     1    -1      0.0000      0.0000      0.0000      0.0000      1.0548     -6.4518      0.0000      1.0600      0.9400
    14     1    -1      0.0000      0.0000      0.0000      0.0000      1.0458     -7.7792      0.0000      1.0600      0.9400
    15     1     1    320.0000    153.0000      0.0000      0.0000      1.0365     -8.0503      0.0000      1.0600      0.9400
    16     1     1    329.4000     32.3000      0.0000      0.0000      1.0468     -6.6353      0.0000      1.0600      0.9400
    17     1    -1      0.0000      0.0000      0.0000      0.0000      1.0466     -7.1527      0.0000      1.0600      0.9400
    18     1     1    158.0000     30.0000      0.0000      0.0000      1.0435     -7.8431      0.0000      1.0600      0.9400
    19     1    -1      0.0000      0.0000      0.0000      0.0000      1.0600     -2.2379      0.0000      1.0600      0.9400
    20     1     1    680.0000    103.0000      0.0000      0.0000      0.9975     -3.1861      0.0000      1.0600      0.9400
    21     1     1    274.0000    115.0000      0.0000      0.0000      1.0450     -4.6395      0.0000      1.0600      0.9400
    22     1    -1      0.0000      0.0000      0.0000      0.0000      1.0600     -0.6449      0.0000      1.0600      0.9400
    23     1     1    247.5000     84.6000      0.0000      0.0000      1.0529     -0.6745      0.0000      1.0600      0.9400
    24     1     1    308.6000    -92.2000      0.0000      0.0000      1.0515     -6.5996      0.0000      1.0600      0.9400
    25     1     1    224.0000     47.2000      0.0000      0.0000      1.0600     -3.0481      0.0000      1.0600      0.9400
    26     1     1    139.0000     17.0000      0.0000      0.0000      1.0596     -4.1867      0.0000      1.0600      0.9400
    27     1     1    281.0000     75.5000      0.0000      0.0000      1.0476     -6.6820      0.0000      1.0600      0.9400
    28     1     1    206.0000     27.6000      0.0000      0.0000      1.0592      0.1778      0.0000      1.0600      0.9400
    29     1     1    283.5000     26.9000      0.0000      0.0000      1.0600      3.1768      0.0000      1.0600      0.9400
    30     2    -1      0.0000      0.0000      0.0000      0.0000      1.0091     -0.7500      0.0000      1.0600      0.9400
    31     3     1      9.2000      4.6000      0.0000      0.0000      1.0600     -0.0000      0.0000      1.0600      0.9400
    32     2    -1      0.0000      0.0000      0.0000      0.0000      1.0280      0.6006      0.0000      1.0600      0.9400
    33     2    -1      0.0000      0.0000      0.0000      0.0000      1.0047      2.3320      0.0000      1.0600      0.9400
    34     2    -1      0.0000      0.0000      0.0000      0.0000      1.0156      2.5267      0.0000      1.0600      0.9400
    35     2    -1      0.0000      0.0000      0.0000      0.0000      1.0583      3.6093      0.0000      1.0600      0.9400
    36     2    -1      0.0000      0.0000      0.0000      0.0000      1.0600      7.2213      0.0000      1.0600      0.9400
    37     2    -1      0.0000      0.0000      0.0000      0.0000      1.0412      3.7914      0.0000      1.0600      0.9400
    38     2    -1      0.0000      0.0000      0.0000      0.0000      1.0395     10.7158      0.0000      1.0600      0.9400
    39     2     1   1104.0000    250.0000      0.0000      0.0000      1.0438     -9.6361      0.0000      1.0600      0.9400
];

% bus       Pg         Qg     Vsetpoint  base   status  Pmax  Pmin  Qmax  Qmin
gen = [
    30    350.0000    -46.3201      1.0091    100.0000     1     1    350.0000      0.0000   9999.0000  -9999.0000
    31    578.9599    362.7264      1.0600    100.0000     1     1   1145.5500      0.0000   9999.0000  -9999.0000
    32    575.2911    224.4590      1.0280    100.0000     1     1    750.0000      0.0000   9999.0000  -9999.0000
    33    563.1107     94.0041      1.0047    100.0000     1     1    732.0000      0.0000   9999.0000  -9999.0000
    34    562.8078    151.5286      1.0156    100.0000     1     1    608.0000      0.0000   9999.0000  -9999.0000
    35    567.7741    200.1804      1.0583    100.0000     1     1    750.0000      0.0000   9999.0000  -9999.0000
    36    564.7114     56.0963      1.0600    100.0000     1     1    660.0000      0.0000   9999.0000  -9999.0000
    37    553.9986     50.2559      1.0412    100.0000     1     1    640.0000      0.0000   9999.0000  -9999.0000
    38    906.5604     48.4800      1.0395    100.0000     1     1    930.0000      0.0000   9999.0000  -9999.0000
    39    970.9868     66.3819      1.0438    100.0000     1     1   1100.0000      0.0000   9999.0000  -9999.0000
];

% fbus  tbus    r        x        b      ratea    rateb     ratec Vf/Vt-ratio  angle  status
% rateA: long term rating; rateB: short term rating; rateC: emergency rating
branch = [
     1     2      0.0035      0.0411      0.6987   9900.0000      0.0000      0.0000      0.0000      0.0000     1
     1    39      0.0010      0.0250      0.7500   9900.0000      0.0000      0.0000      0.0000      0.0000     1
     2     3      0.0013      0.0151      0.2572   9900.0000      0.0000      0.0000      0.0000      0.0000     1
     2    25      0.0070      0.0086      0.1460   9900.0000      0.0000      0.0000      0.0000      0.0000     1
     3     4      0.0013      0.0213      0.2214   9900.0000      0.0000      0.0000      0.0000      0.0000     1
     3    18      0.0011      0.0133      0.2138   9900.0000      0.0000      0.0000      0.0000      0.0000     1
     4     5      0.0008      0.0128      0.1342   9900.0000      0.0000      0.0000      0.0000      0.0000     1
     4    14      0.0008      0.0129      0.1382   9900.0000      0.0000      0.0000      0.0000      0.0000     1
     5     6      0.0002      0.0026      0.0434   9900.0000      0.0000      0.0000      0.0000      0.0000     1
     5     8      0.0008      0.0112      0.1476   9900.0000      0.0000      0.0000      0.0000      0.0000     1
     6     7      0.0006      0.0092      0.1130   9900.0000      0.0000      0.0000      0.0000      0.0000     1
     6    11      0.0007      0.0082      0.1389   9900.0000      0.0000      0.0000      0.0000      0.0000     1
     7     8      0.0004      0.0046      0.0780   9900.0000      0.0000      0.0000      0.0000      0.0000     1
     8     9      0.0023      0.0363      0.3804   9900.0000      0.0000      0.0000      0.0000      0.0000     1
     9    39      0.0010      0.0250      1.2000   9900.0000      0.0000      0.0000      0.0000      0.0000     1
    10    11      0.0004      0.0043      0.0729   9900.0000      0.0000      0.0000      0.0000      0.0000     1
    10    13      0.0004      0.0043      0.0729   9900.0000      0.0000      0.0000      0.0000      0.0000     1
    13    14      0.0009      0.0101      0.1723   9900.0000      0.0000      0.0000      0.0000      0.0000     1
    14    15      0.0018      0.0217      0.3660   9900.0000      0.0000      0.0000      0.0000      0.0000     1
    15    16      0.0009      0.0094      0.1710   9900.0000      0.0000      0.0000      0.0000      0.0000     1
    16    17      0.0007      0.0089      0.1342   9900.0000      0.0000      0.0000      0.0000      0.0000     1
    16    19      0.0016      0.0195      0.3040   9900.0000      0.0000      0.0000      0.0000      0.0000     1
    16    21      0.0008      0.0135      0.2548   9900.0000      0.0000      0.0000      0.0000      0.0000     1
    16    24      0.0003      0.0059      0.0680   9900.0000      0.0000      0.0000      0.0000      0.0000     1
    17    18      0.0007      0.0082      0.1319   9900.0000      0.0000      0.0000      0.0000      0.0000     1
    17    27      0.0013      0.0173      0.3216   9900.0000      0.0000      0.0000      0.0000      0.0000     1
    21    22      0.0008      0.0140      0.2565   9900.0000      0.0000      0.0000      0.0000      0.0000     1
    22    23      0.0006      0.0096      0.1846   9900.0000      0.0000      0.0000      0.0000      0.0000     1
    23    24      0.0022      0.0350      0.3610   9900.0000      0.0000      0.0000      0.0000      0.0000     1
    25    26      0.0032      0.0323      0.5130   9900.0000      0.0000      0.0000      0.0000      0.0000     1
    26    27      0.0014      0.0147      0.2396   9900.0000      0.0000      0.0000      0.0000      0.0000     1
    26    28      0.0043      0.0474      0.7802   9900.0000      0.0000      0.0000      0.0000      0.0000     1
    26    29      0.0057      0.0625      1.0290   9900.0000      0.0000      0.0000      0.0000      0.0000     1
    28    29      0.0014      0.0151      0.2490   9900.0000      0.0000      0.0000      0.0000      0.0000     1
    12    11      0.0016      0.0435      0.0000   9900.0000      0.0000      0.0000      1.0060      0.0000     1
    12    13      0.0016      0.0435      0.0000   9900.0000      0.0000      0.0000      1.0060      0.0000     1
     6    31      0.0000      0.0250      0.0000   9900.0000      0.0000      0.0000      1.0700      0.0000     1
    10    32      0.0000      0.0200      0.0000   9900.0000      0.0000      0.0000      1.0700      0.0000     1
    19    33      0.0007      0.0142      0.0000   9900.0000      0.0000      0.0000      1.0700      0.0000     1
    20    34      0.0009      0.0180      0.0000   9900.0000      0.0000      0.0000      1.0090      0.0000     1
    22    35      0.0000      0.0143      0.0000   9900.0000      0.0000      0.0000      1.0250      0.0000     1
    23    36      0.0005      0.0272      0.0000   9900.0000      0.0000      0.0000      1.0000      0.0000     1
    25    37      0.0006      0.0232      0.0000   9900.0000      0.0000      0.0000      1.0250      0.0000     1
     2    30      0.0000      0.0181      0.0000   9900.0000      0.0000      0.0000      1.0250      0.0000     1
    29    38      0.0008      0.0156      0.0000   9900.0000      0.0000      0.0000      1.0250      0.0000     1
    19    20      0.0007      0.0138      0.0000   9900.0000      0.0000      0.0000      1.0600      0.0000     1
];

% bus   model  startup  shutdown n_coef   Xd'       H=1/2*Tj  D
% model: 1 -- piecewise linear   2 -- polynomial
genpara = [
    30     2      0.0000      0.0000     3      0.0310     42.0000      0.0000
    31     2      0.0000      0.0000     3      0.0697     30.3000      0.0000
    32     2      0.0000      0.0000     3      0.0531     35.8000      0.0000
    33     2      0.0000      0.0000     3      0.0436     28.6000      0.0000
    34     2      0.0000      0.0000     3      0.1320     26.0000      0.0000
    35     2      0.0000      0.0000     3      0.0500     34.8000      0.0000
    36     2      0.0000      0.0000     3      0.0490     26.4000      0.0000
    37     2      0.0000      0.0000     3      0.0570     24.3000      0.0000
    38     2      0.0000      0.0000     3      0.0570     34.5000      0.0000
    39     2      0.0000      0.0000     3      0.0060    500.0000      0.0000
];

% contingency set
% Type    Loctn1   Loctn2    Occurr      Clear_fault    Clear_load     R           X
NTswitch = [
  1110    26    28      0.0000      0.1000      0.3000  1.0000e-010  1.0000e-010
];

% parameters for transient stability
sbspara = [
        0.01
        110/180*pi
       60.00
];
